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Abstract 

We find approximate expressions x{k,n) and y{k,n) for the real and imag- 
inary parts of the kth zero = Xk + iyu of the Bessel polynomial yn{x)- 
To obtain these closed-form formulas we use the fact that the points of well- 
defined curves in the complex plane are limit points of the zeros of the nor- 
malized Bessel polynomials. Thus, these zeros are first computed numerically 
through an implementation of the electrostatic interpretation formulas and 
then, a fit to the real and imaginary parts as functions of k and n is ob- 
tained. It is shown that the resulting complex number x{k,n) + iy{k,n) is 
0(l/'n,^)-convergent to Zk for fixed k. 
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1 Introduction 



The polynomial solutions of the differential equation 

x^y"{x) + 2(x + l)y\x) - n{n + l)y{x) = (1) 

were studied systematically in [T] by the first time. They are named Bessel 
polynomials and are given explicitly by 

where n = 0, 1, ■ ■ ■ . Many properties as well as applications are associated 
to these polynomials: the solution of the wave equation in spherical coordi- 
nates, network and filter design, isotropic turbulence fields, and more (see 
the monograph [2] or [3]-[Tl] and references therein for some other results). 
Among these, several results about the important problem concerning the 
location of its zeros have been obtained [8]-[Il]. Just for instance, the use of 
the reverse Bessel polynomials 6n{x) = x'^yn{^/x) in filter design is known 
since the time when these polynomials began to be studied [21 [H]. Here, the 
poles of the transfer function are essentially the zeros of 9n{x). Thus, it is 
desirable to acquire new analytical knowledge about the location of the zeros 
of the Bessel polynomials. 

In this note we give approximate explicit formulas for both the real and 
imaginary parts of the kth zero Zk = Xk + iyu of yn{x) and show that the 
approximation order of these new formulas to the exact zeros of the Bessel 
polynomials is 0(l/n^) for fixed k. 

The approach followed in this note is simple and based on three items. The 
first is the electrostatic interpretation of the zeros of polynomials satisfying 
second order differential equations [I5]-[I7], the second is a simple curve fit- 
ting of numerical data and the third is the known fact that the points of 
well-defined curves in the complex plane are limit points of the zeros of the 
normalized Bessel polynomials |H]-[11]- The formulas yielded by the elec- 
trostatic interpretation of the zeros of Bessel polynomials are used to find 
them numerically as it has been done previously with these and other sets of 
points [7]-[T9]. Several sets of zeros are computed in this way and the sets of 
real and imaginary values are fitted by polynomials depending on the index 
k whose coefficients depend on n. Finally, it is found that the approximate 
expression for the fcth zero of ?/„(x) is 0(l/n^)-convergent to a limit point 
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of the zeros of the Bessel polynomials. Since the exact zero zj. is 0{l/n'^)- 
convergent to its limit point [S], we conclude that the approximation order 
of our approximate expression to the exact fcth zero is also 0{l/n^). 



2 Asymptotic expressions for the zeros 

Let = Xk + iyki k = 1,2, ■ ■ ■ ,n, be the zeros of the Bessel polynomial 
Unix)- Then, from ([T]) follows that 

N 



k=l ^ 



Zk 



+ 



+ 1 



0, 



z 



where j = 1, 2, ■ ■ ■ , n, i.e., the real and imaginary parts of the zeros should 
satisfy the electrostatic equations 

N 



k=l 

N 

E 



+ 



{xj - Xkf + {Vj - Vkf 



'xj -XkY + ivj -VkY 



+ 



0, 



0. 



(3) 



This set of nonlinear equations can be solved by standard methods. We have 
used a Newton method to solve them up to n = 500. As it is shown in Fig. 1, 
the piecewise linear interpolation of the real and imaginary parts of the zeros 
of the normalized Bessel polynomials yn{x/n) can be fitted by polynomials 
of the second and third degree in the index k. 
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Figure 1: Real and imaginary parts of the zeros of the normalized Bessel 
polynomials ?/„(x/n) for n = 100, 200, 300, 400, and 500, plotted in gray-level 
intensity, from lower to higher, according to the value of n. 
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Thus, we propose the following expressions 

x{k,n) = a2{n)k'^ + ai{n)k + ao(n), (4) 
y{k,n) = b3{n)k^ + b2{n)k'^ + bi{n)k + bo{n) 

to fit our data. To find the relationship between the coefficients of these 
polynomials and n, we take into account the numerical behavior of the data 
at the middle and end points. According to this, x{k,n) and y{k,n) can be 
determined by 

n 3 

x(0,n)=0, ^(2'^) = -^' S:{n + 1)=0, (5) 

and 

~n A 1 + ^ \ n dy{k,n) 4 1 

y{l,n) = , y— — ,n =0, — „ = — , y{n,n) = -. 6 

n 2 ak k=^ n 

These conditions lead to the following coefficients 

6 6(n + 1) 

b (n) - 8(n-2) ^ _ 12(n - 2)(n + 1) 

^2 ^^3 _ 3^2 _|_ 2) ' (n — (n^ — 2n — 2) ' 

, , , 2 (n^ + Qn^ - 12n - A) , , , ra^ - Sn^ + 6 

^iW = -7 .X 0/ 9 — ^ ^oW - 



{n — (n^ — 2n — 2) ' n {n^ — + 2) 

The substitution of ([7]) in (j4]) yields approximate closed-form expressions 

Zk = S:{k,n) +iy{k,n), (8) 

= 1, 2, ■ ■ ■ , n, that converge to the zeros z^. of the Bessel polynomial yn{x), 
as we will show in the following. 



3 Convergence 

Following [9], we define 



and denote by T the curve defined by 



T = {z eC:\W{z)\ = l and |argz| > 

which contains the hmit points Uk of the zeros of the normahzed Bessel 
polynomial yn{x/n). Then, it has been proved in [8] that the zero Uk of 
yn{x/n) approaches to order 0{l/n) the limit value Uk, i.e, 

\ujk-ujk\=0{l/n), (10) 

as n — )■ oo. 

Thus, if we show that |d)fc — (^k\ = 0{l/n), we will have proved that 

- Uk\ = 0{l/n), (11) 

and therefore, taking into account that Uk = nzk, the explicit expression ([S]) 
approaches to order 0(l/n^) the zero Zk of the Bessel polynomial yn{x). 
To this purpose, we simply substitute Cuk = nzk in IQ to obtain, after a 
lengthily calculation, that the expansion of W{uk) in terms of 1/n is 



W{uk) = 1-2 



sin 



Qk 



VlOF -2k + l cos^ f — ^ ] + k 



i + o 

n \n 



,3/2 



4A; -4 
^4:k - 4 

and this implies that 

\Wiujk)\ = l + o(- 

for fixed k. Thus, a)^ approaches to order 0{l/n) the F curve and flTTj) 
follows. From here we have that 

=0(l/n2) (12) 

as n — )■ oo. Numerical calculations confirm and extend this result. Figure 2 
shows the behavior of the maxima of — 5^1 over k as they depend on n. 
The numbers computed by ([3]) are taken as the exact zeros Zk- The displayed 
data shows numerical uniform convergence on the values of k, not only for 
fixed k. A fit of these data gives l/n° with a = 1.7. 
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4 Some few tests 

Just to give examples of the apphcation of the approximate expression ([8]), 
we consider the following cases. 



4.1 The real zero 

An closed-form formula for the unique real zero a in) of the Bessel polynomial 
can be obtained by the substitution of A; = (n + l)/2 in the real part 
of ([8]), x{k,n). This gives 

as our new result. In (21 [TT] asymptotic expressions for a{n) are given. 
Particularly, the following formula 

a(n) = + 0(l/n), (14) 

^ ' 1.32549n + 0.662743 ^ ' ^' ^ ^ 

can be obtained from a more general expression given in [11]. Deleting the 
0-terms and expanding both results in powers of 1/n we find that 

3 1 

a(n) ~ , a(n) ~ -1.50888-, 

2n n 

indicating good agreement between the two approaches. 
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4.2 Power sums 

Here we carry out the corresponding multiplications and use some cases of 
Faulhaber's formula. Then we compare our results with the exact ones. 

1. Sum of the zeros. The simple sum of Zk [cf. ([8])], gives 

^i H = V ^fc = = -1 + O l/n . 

^ — ^ n 

k=l 

The exact result is Si(n) = —1, as can be seen from (|2]). 

2. Sum of the squares of the zeros. In this case we obtain 



S2[n) = y 



E 



q2{n) 
where 



~2 _ P2[n) 

k=l 



p^{n) = 55n® + 15n^ - 800n^ - 612n^ + 4064n^ + 1740n^ 

- 1696n2 - 2832n - 3312, 
q^(n) = 105n^n^ - 2n - 2^ {n^ + n - 2) . 

The exact sum ^ 

^ ' 2n-l' 

can be found elsewhere [12]. Expanding both expressions in powers of 1/n 
we find that 

hin) = ^ + 0(lM, s,{n) = ^ + 0{l/n'). 
21n 2n 

3. Sum of the cubes of the zeros. In a similar form, we obtain that 
where 

P3(n) = -n^° + 28n^-217n*^ + 468n^ + 1002n^- 3804n^- 1076n^ + 5936n^ 

- 3848n2 + 12816n + 19584, 
g3(n) = 35(r^-l)^^^(^^^-6^^-4)^ 

In this case, the exact sum is zero [12j. Expanding ss^n) in 1/n we find that 

h{n) = 0(1/^2). 



7 



5 Final comment 



Note that the approximate formula for Zk given above is not unique. There 
exist other functions to fit the zeros obtained through the electrostatic equa- 
tions (jS]), and there are other conditions to impose at the extreme and middle 
points of the fitting interval. For instance, the imaginary part y{k,n), can 
be fitted by a polynomial of degree 5, but this does not improve the rate of 
convergence and, on the other hand, the calculations become more lengthy. 
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